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ABSTRACT 

In the context of the brane-world scenarios with compactified large extra 
dimensions, we study the production of the possible massive brane oscilla- 
tions (branons) in electron-positron colliders. We compute their contribution 
to the electroweak gauge bosons decay width and to the single-photon and 
single-Z processes. With LEP-I results and assuming non observation at 
LEP-II we present exclusion plots for the brane tension r = and the bra- 
non mass M. Prospects for the next generation of electron-positron colliders 
are also considered. 



PACS: 11.25Mj, ll.lOLm, ll.lSEx 



1 Introduction 



In the last years a lot of attention has been paid to the so called brane world or 
ADD scenarios [1] where the Standard Model particles are confined to live in 
the world brane and only gravitons are free to move along the D dimensional 
bulk space (see [2] for recent reviews). The phenomenological consequences 
of these real or virtual gravitons, described in terms of their corresponding 
Kaluza-Klein (KK) towers, in colliders or in astrophysical and cosmological 
scenarios, have been the object of many recent works (see [3] and references 
therein). However, in addition to the gravitons and the Standard Model 
(SM) particles, one has to consider in principle the possibility of having also 
brane oscillations (branons). In fact, it has been shown that these branons 
or, equivalently, the brane recoil, give rise to an exponential suppression of 
the couplings of the SM particles and the higher KK modes [4] in such a way 
that, in the / <^ Mp regime (where r = is the brane tension and Mp is 
the D dimensional fundamental scale of gravity), the most important modes 
at low energies are the SM particles and the branons. Moreover branons can 
play a role in the solution of some problems appearing when the flexibility of 
the brane is not taken into account such as divergent virtual contributions 
from the KK tower or non-unitary graviton production cross-sections. 

The efi^ective action for the SM fields on the brane was obtained in [5]. 
The introduction of brane fluctuations was done in [6] , where the interaction 
between branons and the SM particles, and the branons self interactions, 
were obtained including also the possibility of having non-vanishing branon 
masses due to the non-factorization of the extra dimension space. 

The effects of brane recoil on real graviton production have been studied 
for example in [7], and on virtual gravitons and gauge bosons in [8, 9]. Some 
constraints from astrophysics on the brane tension were considered in [10] 
and the direct production of branons in colliders was discussed in [11], for 
massless branons in both cases. 

In this work we are interested in the production rates of massive branons 
in electron-positron colliders, in order to get bounds on the brane tension and 
the branon mass in the above mentioned scenario where the brane tension 
scale / is much smaller that the D dimensional fundamental scale of gravity 
Mp. The paper is organized as follows: in Sec. II we define our set up 
and give the effective action for massive branons. In Sec. Ill we obtain 
their couplings to the SM particles and in Sec. IV we discuss the kind of 
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bounds which is possible to set on the brane parameters, namely the number 
of branons n, their mass M and the brane tension scale /. Sec. V is devoted 
to the analysis of the Z invisible width and Sec. VI to the W decay. Direct 
searches based on single-photon and single-Z processes are considered in Sec. 
VII, and extended to future linear colliders in Sec VIII. Sec. IX offers the 
summary and the main conclusions of this work. In addition, the relevant 
Feynman rules are shown in App. A. In App. B the probability amplitudes 
are calculated for the relevant processes and in App. C, it is possible to find 
some two-body phase-space exact integrals which were used in this work. 
Finally, App. D contains some explicit expressions for the cross sections. 

2 The effective action for massive brane fluc- 
tuations 

We will consider a single-brane model in large extra dimensions. In such 
model, our four-dimensional space-time M4 is embedded in a D-dimensional 
bulk space which, for simplicity, we will assume to be of the form Md = 
M4 X B. The B space is a given N-dimensional compact manifold, so that 
D = 4 + N. The brane lies along M4 and we neglect its contribution to the 
bulk gravitational field. The coordinates parametrizing the points in Md will 
be denoted by {x^,y"^), where the different indices run as = 0, 1,2,3 and 
m — 1,2, ...,N. The bulk space Md is endowed with a metric tensor which 
we will denote by Gmn, with signature — , — ). For simplicity, we 

will consider the following ansatz: 



The position of the brane in the bulk can be parametrized as Y = 
{x'^, Y"^{x)), with M = 0, . . . , 3 -|- A*" and where we have chosen the bulk 
coordinates so that the first four are identified with the space-time brane 
coordinates x^. We assume the brane to be created at a certain point in B, 
i.e. Y"^{x) = Y^ which corresponds to its ground state. 

The induced metric on the brane in such state is given by the four- 
dimensional components of the bulk space metric, i.e. g^^, = g^i, = G^y. 
However, when brane excitations are present, the induced metric is given by 



g^.. = d^Y^'d^Y'^GMNix, Y{x)) = ~g^,{x, Y{x)) - d^Y^d,Y-~g'^^{Y{x)). (2) 




(1) 
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Since the mechanism responsible for the creation of the brane is in prin- 
ciple unknown, we will assume that the brane dynamics can be described by 
an effective action. Thus, we will consider the most general expression which 
is invariant under reparametrizations of the brane coordinates. Following the 
philosophy of the effective Lagrangians technique, we will also organize the 
action as a scries in the number of the derivatives of the indiiccd metric over 
a dimensional constant, which can be identified with the brane tension scale 
/. Prom this point of view, to the lowest order in derivatives we find: 

Sb^ [ d'x^ (-f + ...), (3) 

where d^x^ is the volume element of the brane. Notice that this lowest 
order term is the usual Dirac-Nambu-Goto action. However, in certain cir- 
cumstances the effect of the higher order terms must be taken into account 
[12]. 

In the absence of the 3-branc, the metric (1) possesses an isometry group 
which we will assume to be of the form G(M£)) = G(M4) x G{B). If the 
brane ground state is y™(x) = Fq™, the presence of the brane will break 
spontaneously all the B isometries, except those that leave the point Yq 
unchanged. In other words the group G{B) is spontaneously broken down to 
H{Yo), where H{Yq) denotes the isotropy group (or little group) of the point 
Fq- Therefore, we can introduce the coset space K = G{Mjj)/{G{M4) x 
H{Yo))=G{B)/HiYo). 

Let Hi be the H(Yq) generators {i — 1,2,... h), {a — 1,2,... k — 
dim G{B) — dim H) the broken generators, and T = {H, X) the complete 
set of generators of G{B). A similar separation can be done with the Killing 
fields. We will denote those associated to the generators, those 
corresponding to Xq, and by iaiv) the complete set of Killing vectors on B. 
As shown in [6] , the excitations of the brane along the directions of the broken 
generators of B correspond to the zero modes and they are parametrized by 
the Goldstone bosons (branon fields) tt°'{x), which can be understood as 
coordinates on the coset manifold K. When the B space is homogeneous 
then the coset K is isomorphic to B and the isometries are just translations. 
In this case the branon fields (tt) can be identified with properly chosen 
coordinates in the extra space (Y). 

According to the previous discussion, we can write the induced metric on 
the brane in terms of branon fields as: 

QY"^ dY^ 

g,.u = - g'^^^—--jdi,7r''d^7r^. (4) 
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Introducing the tensor hapiT^) as 



hafsin) = /^^;,(r(7r))-^-^, (5) 



we have 



gnu = Qnuix) - -^/ia/3(7r)a^7r"a^7r'^. (6) 

Branons are massless only if the isometry pattern introduced before is ex- 
act. However, in general, the symmetry is only approximately realized and 
branons will acquire mass. In order to show this mechanism explicitly, let us 
perturb the four- dimensional components of the background metric and let 
g^i, depend, not only on the x coordinates, but also on the y ones: 

r< I 9f.u{x,y) \ , . 

- [ --g'^^iy) ) ■ 

This has to be done in such a way that the G{B) piece of the full isometry 
group is explicitly broken. Notice that the breaking of the G{B) group by 
perturbing only the internal metric g'mniv) '^^^^ not lead to a mass term for 
the branons. 

In order to calculate the branon mass matrix, we need to know first the 
ground state around which the brane is fluctuating. With that purpose, we 
will consider for simplicity the lowest-order action, given by: 

^S/N = / d'xC^'^ = -f[ d'x^g{x,Y{x)) (8) 

which will have an extremum provided 

5sfj^[K] = ^ = = ^ r^'dm-g,. = 0, Vy-. (9) 

This is a set of equations whose solution Y^{x) determines the shape of the 
brane in its ground state for a given background metric g^i,. In addition, the 
condition for the energy to be minimum requires: 



52£(0) 



Y=Yo 



f^f^r^idndmg^cu - '2r''dn~g,admgi.p) < 0. (10) 



i.e. the eigenvalues of the above matrix should be negative. This implies 
that the static Lagrangian should have a maximum. 
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In order to obtain the explicit expression of the branon mass matrix, we 
expand g^u{x,y) around = FJ" terms of the tt" fields: 

~9A^,y) = ~g^,.{x,Yo) + dm~9^.,{x,Yo){Y"'-Y,"') (11) 
+ ldmdn~g,uix,Yo)iY^ - Y,^)iY^ - V) + ... 

The linear term in branon fields is written as: 



1 f'" 

jV^l}u-dm~g,.{x,y)\^^y^p, (12) 



while the quadratic term takes the general form: 



TT = 



(13) 



(3) 

and we have not written the explicit expression for V^y^^j, since it will play 
no role in the present work. Here, we have used the fact that the action of an 
element of G{B) on B will map Yq into some other point with coordinates: 

Y^ix) = Y^iYo, n^ix)) = FJ" + ^^^(Yo)7r-{x) + 0{n') (14) 

where we have set the normalization of the branon and Killing fields with 
= IGn/Mp being Mp the four-dimensional Planck mass. Substituting the 
above expression back in (6), we get the expansion of the induced metric in 
branon fields: 

g^,^~g^,(x,Yo)-j^5o.pd,n-d.n^ + jV^l},n'^ + j^vi^^^^^ (15) 

We have also used the fact that since vr" must be properly normalized 
scalar fields, the Y"^ coordinates should be normal and geodesic in a neigh- 
borhood of Y^ and, in particular, they cannot be angular coordinates. This 
implies that we can write ha/si'K = 0) = Sa/s- 

Assuming for concrctcness that, in the ground state, the four-dimensional 
background metric is flat, i.e. ^^j^(,x, "Ko) = Vfn^^ the appearance of the 
^aia2---aifj.iy tsusors in (11) could break Lorentz invariance, unless they fac- 
tor out as = M^}a2...aiVfiu/{4:p)- With this assumption, the linear 
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term vanishes identically due to the condition of minimum for the brane 

energy (9), and the M^^ coefficient in the quadratic term can be identified 
with the branon mass matrix. Thus, in general and for the square root of 
the metric determinant we find: 

= 1 - ^^'^'^a/ja^TT^a.TT'' + ^y^y^^ + - (16) 

Notice that this expression requires that both dn / and M^n^/f^ arc small. 
This includes different types of expansions, such as low-energy expansions 
with small branon masses compared to /, or low-energy expansions with 
possible large masses and small n/f factors. 

The different terms in the effective action can be organized according to 
the number of branon fields: 

Seff[n]^Si%[n]+Si%[n] + ... (17) 

where the zeroth order term is just a constant. The free action contains the 
terms with two branons: 

•^S/M = U d^4^^f,d,n'^d^n^ - Ml^n'^n^). (18) 

Z, J Ma 



3 Couplings to the Standard Model fields 

As we have shown in the previous sections, the induced metric on the brane 
depends on both the four- dimensional bulk metric components and the 
branon fields vr". As noticed in [8], in the limit in which gravity decouples 
Mp oo, the branon fields still survive. This implies that their effects can 
be studied independently of gravity. With that purpose, in the following, we 
will consider their couphngs to the Standard Model fields in the absence of 
gravitational background field i.e. g^^i, = rj^^. In order to obtain the general 
couplings, one can proceed as in [6], where the action on the brane, which is 
basically the SM action defined on a curved space-time M4, is expanded in 
branon fields through the induced metric. For example, the complete action, 
including terms with two branons, contains the SM term, the kinetic term 
for the branons and the interaction terms between the SM particles and the 
branons: 



6 



Sb ^ [ d^'x^i-f + CsMig,.)] 

J Mi 



j ■ 



0{^% (19) 



where Tg'^^rj^iy) is the conserved energy- momentum tensor of the Standard 
Model evaluated in the background metric: 



rplJ,V 



5C 



SM 



(20) 



It is interesting to note that there is no single branon interactions which, 
as commented above would be related to Lorentz invariance breaking. In 
addition the quadratic expression in (19) is valid for any internal B space, 
regardless of the particular form of the metric In fact, the form of the 
couplings only depends on the number of branon fields, their mass and the 
branc tension. The dependence on the geometry of the extra dimensions will 
appear only at higher orders. 

Let us give the results for the different kinds of fields. The corresponding 
Feynman rules can be found in Appendix A (we will follow the notation in 

[6]): 



Scalars 

For a complex scalar field $ with mass m$ in a certain representation of a 
gauge group with generators T", the Lagrangian is given by: 

Cs = (D^^YD^^ - m|$t$, (21) 

whose energy-momentum tensor is given by: 

T^" = -ri'"'{{Dp^)^D''^ - m|$t$) + {D^'^)^D''^ + {D'^fSy^D"^ (22) 

where the gauge covariant derivative has the usual form = — /lA^T". 
From this expression we find the Feynman rules for the following types of 
vertices: 7r7r$^$, nn^^^A and nn^^^AA. 
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Fermions 

The Lagrangian for a Dirac fermion field with mass reads: 

Cf = il;(trD^-m^)tlj, (23) 
and its associated energy-momentum tensor is: 

- \n''''{i{i>YDp-Dpi>Y)i^-'^rn^m- (24) 

where again we assume the fermion field ip to be in a certain representation of 
a gauge group with generators T", and the covariant derivatives are defined 
as = dn — hT"-A1^{cv{a) —0^(0)75), where, in general, the vector and axial 
couplings could be different. The Feynman rules will contain the following 
types of vertices: Tnr'ipip and inr'ipipA. 

Gauge bosons 

For the Yang-Mills action on the brane we can follow similar steps. The 
Lagrangian is given by: 

jC^^Af-^-F;, + Up (25) 

where: 

F;, = d^Al - d^Al - hC'^^A'^Al, (26) 

jCfp is the Fadeev- Popov Lagrangian, which includes the gauge fixing and 
the ghost terms. The gauge symmetry can be spontaneously broken in such 
a way that the gauge fields acquire masses through their couplings to 
the scalar sector (21). In that renormalizable gauge will be used in 

the calculations. It is interesting to notice that the metric used in the gauge 
fixing term could be either the induced metric or the flat background metric. 
Both choices give rise to valid gauge-fixing expressions. The second one is 
however simpler because it does not introduce new couplings between gauge 
or Goldstone bosons fields and branons. The same criterion has to be taken 
for the ghost Lagrangian, which will not be studied here since we are only 
interested in the tree-level analysis. 
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The energy-momentum tensor takes the form: 

T7 = F;,F«,(r?'^VV^ + ^^'VV^)+T^P (27) 

and from this interaction action, one finds the Feynman rules for vertices hke 
-K-nAA, nnAAA and -n-nAAAA 

We see that in all these vertices, branons always interact by pairs with 
the SM matter fields. In addition, due to their geometric origin, those in- 
teractions are very similar to the gravitational ones since the tt fields couple 
to all the matter fields through the energy-momentum tensor and with the 
same strength, which is suppressed by a factor. In fact, branons couple 
as gravitons do, with the following identification: 

V -74/ia/3(7r).9^7r"5,7r^ (28) 

where /i^,^ is the graviton field in linearized gravity. 



4 Constraining brane-world models 

In the brane-world scenario, with / ^ Mp^ the relevant low-energy excita- 
tions on the brane correspond to the branons rather than to the Kaluza-Klein 
graviton excitations. In such case, and provided the energy scale of present 
accelerators is well below both the fundamental scale of gravity Mp and the 
brane tension scale /, but is large compared to the branon mass M, such 
branon fields will be produced in the collision of SM particles. Moreover, 
the calculation of the corresponding production cross-sections could be per- 
formed by using the effective theory described in the previous section. 

In principle, there would be new physics related to both the KK gravitons 
and the branons, but the dependence on / in each case is completely differ- 
ent. For example, the behavior of typical cross sections for missing energy 
processes, i.e. KK graviton and branon production processes, becomes large 
for both small and large values of the brane tension / [10]. The cross section 
for a small value of / is governed by the branon production process, while for 
a large value, the KK gravitons production rate dominates. Therefore if the 
brane tension scale is smaller than the fundamental gravitational one, then 
the first indications of extra dimensions would be given by the production of 
branons. This would allow us to measure the brane tension and the branon 
mass or, if branons are not found, at least to set bounds on such parameters. 
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The effective couplings introduced in the previous section provide the 
necessary tools to compute the cross sections and the expected rates of new 
exotic events in terms of / and the branons masses only. In fact, with a 
rotation in the coset space K, the branon mass matrix takes a diagonal 
form. We will assume that there are n branons with the same mass M. If 
the branons masses are similar this would be a good approximation. In the 
opposite case, with very different masses, we could neglect the heavier fields 
and consider only the production of light branons. In these simple cases, all 
the cross sections will be parametrized by three unknown parameters: 

n: number of branons. 
M: mass of the branons. 
/: brane tension scale. 

In the following, we will compare the data coming from LEP with the pre- 
dicted cross section for different processes in which we expect high sensitivity 
to new physics. In particular we will concentrate on the branon contribution 
to the invisible width of the Z boson, to the decay rate of the bosons and 
to direct searches with single-photon and single-Z vertices. This will allow 
us to set bounds on certain combinations of the above parameters. 

5 Bounds from the Z invisible width 

In the Standard Model, the full decay width of the Z boson has three types 
of contributions, coming from charged leptons, hadrons and neutrinos re- 
spectively. The third one corresponds to the so called invisible width, since 
neutrinos escape without producing any signal in the detectors. In princi- 
ple, physics beyond the SM could also give rise to invisible decays of Z and 
therefore to deviations with respect to the SM predictions. 

The precision measurements done mostly at LEP-I set stringent limits 
on the Z invisible decay width, which can be translated into strong bounds 
on the new channels contributing to such process. In fact in the following, 
we will show how it is possible to get lower bounds on the brane tension / 
coming from light branons (M < in which the Z could decay. Since 

branons are stable and chargeless particles, they would also avoid detection. 

In this section we calculate the first order contribution from branons 
to this width, which is given by the decay of Z into two branons and two 
neutrinos: 

r|: Z — ^ t>(pi),t^(p2)7r(A;i),7r(/c2) 
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Wc use the expression for the probabihty amphtudes in Appendix B 
Eq.(63) calculated with the Feynman rules given in Appendix A. We in- 
tegrate over the branon two-body phase space making use of the formulae in 
Appendix C. The result is averaged over the three Z polarizations. In addi- 
tion, such expression should be divided by two because the outgoing branons 
are indistinguishable. Thus we get: 

dr^ _ 2Mzn 



+ 
+ 

+ 
+ 

+ 
+ 

where we have used equations (63) with cy = = 1, Ma = Mz and the 
coupling constant h = g/{4:Cos9w), with 9w the Weinberg angle. We have 
also made use of the Mandelstam variables: s = (pi + ]?2)^, t = {pi — qf, 
u = {p2—qY and k'^ = (A;i+A;2)^, although only three of them are independent 
because s + t + u = k'^ + M|. 

We can also compute the above differential cross section in terms of (non- 
covariant) variables with a more physical interpretation, namely, the energies 
of the outgoing neutrino and antineutrino, Pi = Ei and P2 = E2 respectively, 
and the angle between their three- momenta 9. In terms of these variables we 
have: 

s = 2p?p0(l-cos^), 
t = Mz(Mz-2p°), 

u = MziMz-2pl), (30) 
The total width is calculated by integrating the above expression with respect 
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W 1 - -p-{20M^M|(2A;^ - 5M2)tV(M|(2s - A;^) + tu) 

(A;2 - 4:M'^ f{stu{s{k'^ + M|) + Uu){2s{k^ + M|) + u^) 
{t^ + u^){2s + 2e + Ml)Ml 

[6s{t^ + u^) + 6s^{t^ + u^- tu) + 3tu{t - uf + t^ + u^]Ml 
M|[2s^(2(t2 + u'^) - 5tu) + 2s2(3(t^ + u^) - 5tu{t + u)) 
s(2(i^ + u^) + tu{t^ + - 8tu)) + tu{t^ + - 7tu{t + u))] 
Ml[s\t -uf- 8iV + 2s^{t^ + u^- 2tu{t + u)) 
2stu{t^ + u^- Stu{t + u)) 

s^t^ + u^ + tuit^ + u^-Utu))]}} (29) 



to the three variables within the following kinematic limits: 

v\ e 7,.[0,^^], (31) 



^ Mz{Mz-2p\)-AM^ 
' 2(Mz- (l-cos%?)^' 
cos^ e 7o=[-i,i]- 



The final result, depending on the three undetermined parameters reads: 

r|(n,/,M)=/ „f dp\dpld cose. (32) 

The integration over the rest of angles has been done immediately because 
the decay rate does not depend of them. It is interesting to note that the 
dependence on the number of branons n and brane tension / factorizes, and 
only the dependence on M is non-trivial: 

n 

r|(n, f,M) = ^x U{M/Mz) (33) 

where n(a;) is a dimensionless function, which depends on a single dimen- 
sionless variable, and which can be easily calculated from (32). For example, 
for massless branons a numerical integration gives: 

r'z(^,/,0) = ^x6.15xl0-^^ (34) 

To show the dependence on the branon mass, we plot in the Figure 1, 
rz(^, /, ^)/^/^ against M. Finally, to get the total Z invisible width into 
branons, we have to multiply by a factor of 3 because there are three different 
neutrino families. 

The key observation is that at the 95% confidence level, the variation of 
the Z invisible width cannot be larger than 2.0 MeV [13], i.e., the contribution 
from new physics satifies: 

r5^'"(n,/,M) < 2.0 MeV, (35) 

It is then possible to find bounds for the different brane parameters just 
imposing the above limit on Eq.(33). For example, if branons are massless, 
the bound on / depends only on the number of branons: 

/ > ll.gni/^GeV. (36) 

On the other hand, in the case of one branon, we show the exclusion plot 
in the f — M plane (see Figure 2). 
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M (GeV) 

Figure 1: and Z widths as a function of the branon mass. Both plots correspond to 
a single channel. We have extracted the dependence on the brane tension and the number 
of branons in the factor /n. 

6 Bounds from decay 

There are mainly two channels contributing to this process. First, the decay 
of into two branons and two leptons: 

T^: W- ^l-ipi),iyip2),n{h),n{h) 

Such leptons can be an electron and an electron antineutrino for the W~ 
decay, their antiparticles for W^, or the analogue pairs of leptons in the rest 
of families. The results of the different decay rates agree in the limit in which 
the lepton masses are negligible. In fact this is a good approximation since 
the typical energy carried by the of outgoing particles will be comparable to 
the W"^ mass which is much larger than the leptons masses. 

Second, we have the W"^ decay into two branons and a quark- ant iquark 
pair: 

W — > q{pi), q{p2), rc{ki), 77(^2) 

In principle there should be also a channel including gluons in the final state, 
however in the SM, such channel is suppressed by a coefficient as{Mw)/7r = 
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I Unexcluded area by the invisible Z width at the 95% CL. 

^1 Excluded area by the invisible Z width at the 95% CL. 
^] Excluded area by the visible W width at the 95% CL 

Figure 2: Exclusion plot in the f — M plane for single-branon models from LEP-I data. 
The dark area is excluded by the measurement of the Z invisible width and the striped 
region is excluded by the measurement of the visible decay width. 

0.04 and for that reason we will also ignore it in the present case. In any 
case, the result we will obtain will be a strict lower bound to the true decay 
width into branons. 

The calculation is totally similar to that of the Z decay and for that 
reason we will not repeat it here. The expressions for the Z decay in (29) 
are also valid for the W^, the only changes, apart from the the different 
gauge-boson mass, are the couplings to leptons and quarks. Thus, in 
this case we have: Ma = Mw, cy = ca = I, h = g/{2\/2) for leptons and 
h = Vijg /{2\/2) for quarks where is the Cabbibo-Kobayashi-Maskawa 
(CKM) matrix element with i = u,c,t and j = d, s, b. The differential decay 
rate has the same form as in the Z case, which is written again in terms of 
the energy of the two leptons or the quark- antiquark pair p^, P2 and the angle 
between their three momenta that we denote by 6. Thus the corresponding 
Mandelstam variables are obtained from (30) just replacing Mz — >■ Mw The 
total integration over the leptons or quarks phase-space has to be done in a 
numerical form, with the analogous limits to those in (32) again replacing 
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the gauge bosons masses. 

The final result r^(n, /, M) depends on the three unknown parameters 
and is obtained from the differential rate by using the analogue expression 
to (32). We see that once again, the dependence on n and / factorizes, and 
the dependence on M enters through the same dimensionless function Il[x). 
Thus we can write for the leptonic decay: 

C(n,/,M) = 2cos^9w^ x U{M/Mw). (37) 

For massless branons we obtain by numerical integration: 

■n 

^win, f, 0) = = X 9.56 x 10"^^ (38) 

The results in the case of quarks can be obtained in a straightforward way 
from the leptonic one just multiplying by the modulus squared of the corre- 
sponding CKM matrix element. 

To show the dependence on the branon mass, we have plotted T^f^/n 
against the branon mass M in Figure 1. 

The most important difference with respect to the Z analysis arises from 
the experimental constraint on the width coming from new physics. In 
this case, the decay is visible and, for example, if the final decay products con- 
tain e~ and Pg, the process with branons could give a signal similar to the fol- 
lowing Standard Model process: F^: W~ — >• e~{pi),Ue{p2),^T{P3),i^T{P4), 
in which, the decays into a r~ and a P^, and then the r~ decays finally 
into an e^, and z/^. 

The behaviour of the W"^ decay width in terms of the visible lepton energy 
El is plotted in Figure 3 for different branon masses. 

However in this work, we are only interested in the constraints on the 
theory parameters coming from the total width. The total leptonic 
decay is obtained multiplying (37) by a factor of 3, because there are three 
different kinds of processes corresponding to the three leptonic families. On 
the other hand, to calculate the hadronic decay we have to multiply by: 

3 J2 I V^ijf = 6.117 ±0.075, (39) 

i=u, c ;j=d, s, b 

where the factor of 3 comes from the numbers of colors and the sum is 
extended to the six CKM matrix elements which do not involve the top 
quark. The numerical result is taken from the combined LEP measurement 
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Fi gure 3: Behaviour of the decay width in terms of the visible lepton energy Ei for 
different branon masses. 

of Br{W — > qq) [14]. The above approximation assumes that the mass 
of these five quarks is neglected when compared to the typical energy of 
the process {Mw), and the top contribution is neglected because of its large 
mass. We can thus estimate the total width just multiplying by 3 + 6.12 
the result in (37). 

In order to get bounds on the contribution from new physics to the 
width, we take into account the uncertainties on the experimental 

measurement of the visible W width: = 2.150 ± 0.091 GeV [13], and 

the uncertainties in the SM calculation: = 2.093 ± 0.003 GeV [14]. So 
that for the total allowed variation in the visible width we obtain at the 
95% confidence level: 

r^V^"'(n, /, M) < (r^^^ - T^J') + 2Ar^^^ = 240 MeV, (40) 

This limit translates into the following one for massless branons: 

/ > 6.9^1/^ GeV. (41) 

The exclusion regions as a function of the variables / and M are shown 
in Figure 2 for the single branon case. 
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7 Branon bounds from direct searches 



7.1 Cross section with a generic gauge boson 



In this section we present the cross section of the process: 



n{h),n{k2),Al{q) 



(42) 



The initial particles can be either leptons or quarks, although we will be 
mainly interested in the e'^e~ case. There are three final particles, one ar- 
bitrary gauge boson (either Z or 7) together with the two branons. We are 
interested in the differential cross section with respect to the gauge boson 
phase-space parameters, i.e. energy and scattering angle. 

We have integrated over the two-body phase space of branons. The prob- 
abilities are calculated as usual by doing the square of the amplitude module 
given in (62). We have also averaged over the spin of the initial particles and 
summed over polarizations of the outgoing boson. The probability has to be 
divided by two because the outgoing branons are indistinguishable. 

We define the invariants: s = (pi + ^2)^, t = {pi — q^jU = {p2 — qY 
and k'^ = {ki + A;2)^- Again only three of them are independent because 
s + t + u — k'^ + {Ma is the gauge boson mass and we are neglecting 
the fermion masses). In the CM system the probability should be divided by 
2p°2p>i - W2I = 2s. 

Again we can change to more physical variables, namely, the energy frac- 
tion of the emitted gauge boson x — 2go/-\/i and the scattering angle 9. Thus 
we get: 



u 



t 



-V^iqo-\t\cosi9)) + M', 
-v^{qo + \t\cos{9)) + M^ 

Xy/S 



(43) 



2 




and the phase-space volume takes the form: 




dk^dtdcf) dxd{cos9)d(i) 



^J{sxf-AsMl. (44) 



2go(27r)3 



4s(27r)3 8(27r)3 
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Therefore, for Ma ^ and arbitrary cy and ca parameters, we have: 



da A 
dxd cos 9 



LU2 



47r 122880/%2(s_M2)2^ 



2 - 4M2){5M2m2(2(s(1 - x) + M^) - 5M2) 



{8M„2 + (sx^ - AMI) sin^ ^} 

2(s(l - x) + Ml - AM^f 



+ 



(4M2(M2 + s(l - x)) + s(sx2 - 4M2) sin^ Of 
{lQMl{s{l -x) + Mlf[lQMl + M„S(9 - 5x) 
+ 2M^s2(4-3x + x2) + s=^(l-x)] 
- 2[40M„^2 - 2Mf s(-112 + 40a; + 5x^) 
+ 4M^s^(78 - 70x - + 5x^) + s^x^(-2 + 4x - + x^) 
+ Mfs^(216- 272x + 6x^ + 66x^-llx^) 
+ M„^s^(96 - 184x + 62x^ + 44x^ - 21a;^ + x^) 

M^s^iS - 16x - I2x^ + A2x^ - 29x^ + 6a;^)] sin^ 6 



+ 
+ 



s{sx' - 4M^y[7M^ + M^s(17 - 7x) + M>'(1 - 2a;) 



+ s%3{l -x) + 2x')] sin^ ^ - s\sx' - 4.M^y{s - M^) sin^ 6}}. 

Notice that this expression is also apphcable to quark-antiquark coUisions 
just by choosing the appropriate couphngs (cy and ca parameters together 
with h). 

Also, for certain processes it is interesting to calculate the cross sections 
without summing over the polarizations (ep of the outgoing gauge boson. 
Again we consider the case Ma ^ and arbitrary cy and ca parameters. The 
results can be found in Appendix D for longitudinal (72) and transverse (73) 
polarizations. It is interesting to note that the longitudinal cross secction 
vanishes in the limit in which the gauge boson mass goes to zero. Therefore, 
the total cross section (including the three polarizations) in the limit M<j — 
has the same form as the total cross section for a massless gauge field 
{Ma — 0), and it reads: 



da. 



\h\^s{4 + c\){s{l-x)-4M' 



2\2 



n 



dxd cos 9 



47r 61440/%2 

a;(3 - 3a; + 2a;^) - a;^ sin^ 9 + 



A 



1 - 



4M2 



s(l — x) 

2(1 -a;)(l + {l-xf) 
X siv? 9 



(45) 
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As expected, the cross section is divergent for collinear outgoing boson [6 = 
0) and also for a soft {x = 0) massless gauge boson. 

7.2 Bounds from single-photon processes 

The differential cross section obtained in the previous section could be used 
in direct searches of branons in colliders. In this section we are going to 
calculate the contribution to the total cross section of processes involving a 
single photon in the final state. This will allow us to obtain new bounds on 
M and / assuming non observation at LEP-11. 

To perform the angular integration over the polar angle cos^ in (45), it 
is necessary to take into account the angular range covered by the detector 
cos^ e [—d,d\. Thus we get: 



dx 



+ 

where h = e, cy = I and ca = 0. 

The integration in the x variable is done within the following limits x e 
[Ejn/ {2^/s), l—AM^/s]. The upper bound is fixed by kinematical constraints, 
but the lower one {X = E„J2^), imposing a minimum energy for the 
photons (Em), depends on specific experimental conditions (trigger, noise, 
backgrounds, . . . ). To perform the numerical study we take d = 0.96 and 
Ejn — 5 GeV, which correspond to the typical boundary conditions of one 
of the LEP experiments (see [15]). We will also use ^/s — 206 GeV to do 
the calculations. Although the maximum centre-of-mass energy was close to 
209 GeV, most of the data at LEP-11 taken in the last year of running were 
collected in the range between 205 and 207 GeV [13]. 

The total cross section can be calculated analytically, although it will not 
be shown here since the final expression is quite long. The result for massless 
branons (M=0) is much simpler, being just proportional to s^: 

<^t(^ = 0) = (47) 
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\h\'s{cl + c\)isil 



x 



AM^Yn 



1 - 



47r 184320 /8 7r2 a; \ 

{2dx^ (9-9x+ (3 + ci^) x^) 

3 (-2 + Ax-?>x^ + x^) log f^^— ^ 



4M2 
sil — x) 



(46) 




M(GeV) 

Figure 4: Total cross sections for single-photon and single-Z processes as a function of 
the branon mass. The dependence on the brane tension and the number of branons has 
been factorized out in the /^/n coefficient. 



where the constant depends on the detection hmits as follows: 



A, 



^ 11059200 7r2 

+ (^1 + 3X + 6X2 + 10X^ 



2 d (-1 + X f [30 ( 1 + 3 X - 3 + 



+ 3 log 



Xl+df 

- 390 + 220 X^ - 75 X^ + 12 X^ - 120 log(X)' 



-247 + 480 X 



(48) 



with aEM = e^/(47r). 

For the limits mentioned before we get A^ = 4.54 x 10~^. The branon 
mass dependence of the total cross section is plotted in Figure 4. 

To obtain the bound on the theory parameters, we can consider the limit 
on the total cross section coming from new physics which could be added 
to the single photon channel without being detected. For LEP-II, this limit 
is approximately anew — 0.1 pb, which is the typical experimental sensitiv- 
ity observed in LEP searches [16]. The result for one branon assuming no 
observation is plotted in Figure 5. Thus the interior area limited by the 
bound curve is potentially excluded by the LEP-II experiment. In the case 



20 



/ (GeV) 




Unexcluded area by the single photon signal. 
Excluded area by the single photon signal. 
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Figure 5: Exclusion regions in the f — M plane coming from single-photon events (dark 
region) and single-Z processes (striped area), using LEP-II data. 



of massless branons, the limit depending on the number of branons reads: 

/ > 138ni/^GeV. (49) 



7.3 Bounds from single-Z processes 

As in the single-photon case, the single-Z channel can be used to restrict the 
parameters of the theory. In this section we are going to calculate the total 
cross section summing over Z polarizations to estimate these bounds. 

Again to implement the angular integration over the polar angle cos^, 
we take into account the domain of the detector: cos^ G [— c?, rf], with d = 
0.96. The behavior of the cross section in terms of the outgoing Z energy is 
represented in Figure 6 for different branon masses (for y/s = 206 GeV). In 
this case, h = g/{4:Cos9w), 9 = sm9w, cy = — 1 + Asin^ 9w and ca = 1- 
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Figure 6: Differential cross section for single-Z processes as a function of the outgoing 
Z energy for different branon masses and for a center of mass energy y^s = 206 GeV. The 
dependence on / and n has been factorized out. 

The integration in the x variable is done within the following limits x G 
[Mz/i^y/s), 1 + (M| — 4M^)/s], where both the lower and the upper limits 
are kinematical. 

To estimate a bound on the theory parameters, we will consider again the 
limit on the total cross section coming from new physics that can be added 
to the single Z channel without being detected. For LEP-II, the order of 
magnitude of this limit is approximately equal to that in the single photon 

case CTfiew 

~ 0.1 pb. 
The result for massless branons is: 



and the exclusion plot can be found in Figure 5. The interior area limited 
by the bound curve, which is excluded, is smaller than in the single photon 
analysis. This is an expected result since the Z coupling to the electron field 
is smaller than that of the photon, and the Z mass restricts the avalaible 
phase space. In fact, this restricts the search only to branons with masses: 
M < (a/s — Mz)l1 — 57 GeV. In the single-photon channel however, the 
kinematical range is larger M < \/s/2 ~ 103 GeV. Despite this fact, the 



/ > 52n 



i/«GeV. 



(50) 
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single Z channel analysis is still interesting since it provides a completely 
independent direct search method. In future colliders, working at very high 
centre-of-mass energy, the Z mass could be neglected ^/s » Mz, and the 
total cross section for the two channels will take the same form. The only 
difference will come from the ratio of couplings to the electron field, which 
is: 

^(V5»A«.il!%^i?±i.0.37, (51) 

a-y 16 cos^ Ow sm Bw 

where we have use the value sin^ 9^ — 0.22. This implies that provided the 
rest of the analysis remains unchanged, the bounds on / which come from 
these two processes are similar: 

^^{V^ » Mz) ~ (0.37)1/8 ^ Q_88^ (52) 

Although, in general, the single-Z channel has more background compared 
to the single photon case, implying lower precision, however it allows us to 
perform analysis depending on the polarization which could improve the 
bounds. 

Since branon effects grow strongly with energy, it is not surprising that the 
bounds from direct searches with {^/s ~ 200 GeV) are more constraining than 
the indirect ones, in which the energy scale is set by the Z mass {Mz ^ 90 
GeV). However, because this is an effective theory the growth with energy 
of the cross sections will eventually violate unitarity. At that point the 
approximation will no longer be valid. It is also interesting to note that the 
present bounds improve the astrophysical ones for massless branons, coming 
from the 1987a supernova, from which / > 10 GeV [10]. 



8 Prospects for future linear colliders 

Several proposals for the construction of e+e" linear colliders in the TeV 
range are currently under study. The TESLA (TcV Energy Superconducting 
Linear Accelerator) [17], the NLC (Next Linear Collider) [18] and the JLC 
(Japanese Linear CoUider) [19] are examples of the first generation of these 
colliders, whereas the CLIC (Compact Linear Collider) [20] would correspond 
to the second generation. In this section we discuss the sensitivity of these 
colliders to a hypothetical branon signal. The study will be performed in 
terms of the brane tension, /, and the branon mass, M. 



23 



The physics programme of the new hnear colhder projects includes the 
measurement of electroweak parameters with improved precision, such as the 
invisible Z width or the width. However, and since the deviations due 
to the presence of branons increase dramatically with energy, the largest 
sensitivity to a branon signal is expected in direct searches like single-7 and 
single- Z. 

In order to estimate the sensitivity of future linear colliders to branon 
signals, the LEP-II study of the previous section is extended to higher center- 
of-mass energies. It is assumed that, at the time of construction of these ac- 
celerators, the theoretical and systematic uncertainties on Standard Model 
processes will be controlled at the level of the femtobarn. Under this as- 
sumption, we will estimate the sensitivity limit, cf^^w i ^"7 scaling the LEP-II 
estimate by the expected gain in statistics: 



^FLC 
new 



{TILleP-Ii) n 1 u (r.o\ 



where TILlep-ii — 700 pb~^ is the LEP-II integrated luminosity and 0.1 pb 
is the LEP-II sensitivity limit. We will consider the following values of the 
integrated luminosity at future colliders, TILpLC'- TILi ~ 200 fb~^ for a 
first stage of TESLA, NLC or JLC, and TIL2 — 1000 fb~^ as a maximum 
value for a second stage. For CLIC we will assume the same total integrated 
luminosity: TILclic = TIL2 — 1000 fb^^ The cross section bounds for 
both luminosity choices are similar: cr^^g^ = 6 fb and cr^g^ = ^ntw^ = 3 fb. 

The critical parameter in the analysis is the center-of-mass energy, Eqm- 
In the single photon channel, this leads to limits on / for any branon mass 
M < Ecm/'2- In particular the limit for massless branons (/q) increases pro- 
portionally to according to Eq. (47) . Assuming a center-of-mass energy 
for the first stage of the first generation of hnear collider of approximately 500 
GeV, 1 TeV for its second stage and 5 TeV for CLIC, we obtain the following 
limits for /o in the case of a single branon: f^^ > 398 GeV, /(f^*^ > 758 GeV 
and /q^^^ > 2.64 TeV. The results of a full study in the (/, M) plane and 
in different experimental contexts are presented in Figure 7. 

For the single-Z channel, the bounds are less restrictive: /q^ > 205 GeV, 
f^z > 450 GeV an f^^^'^ > 1.87 TeV. Obviously, the study in this case is 
only applicable to branon masses below {Ecm — Mz)/2, due to kinematic 
constraints. The excluded regions in the (/, M) plane are also shown in 
Figure 7. 
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Figure 7: Predicted experimental accessible regions in the / — M plane via single-photon 
and single-Z processes. Different center-of-mass energies {Ecu) and total integrated lu- 
minosities (TIL) at future e+e^ linear colliders are considered. 
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9 Summary and conclusions 



In this work a brane- world scenario where the brane tension scale / is much 

smaller than the fundamental gravitational scale Mp has been considered. 
For this case, the relevant low-energy degrees of freedom are the brane SM 
particles and the brane oscillations or branons. From the corresponding ef- 
fective action, and including also the effects of a possible branon non-zero 
mass, we have obtained the relevant Feynman rules for the couplings of bra- 
nons to SM particles. They have allowed us to compute the decay rates and 
cross-sections for the different processes relevant for branon production in 
electron-positron colliders. We have used the information coming from LEP 
in order to get different exclusion plots on the branon mass and the tension 
scale plane. Single-photon production turns out to be the most efficient pro- 
cess in order to set bounds on the brane parameters / and M. We have 
also extended the analysis to future electron-positron colliders. The corre- 
sponding exclusion plots that could be obtained in case that branons were 
not observed, have been also shown. 

The work presented here should be complemented with a parallel analy- 
sis for electron-hadron colliders such as HERA and hadron-hadron colliders 
such as the Tevatron or the LHC, and with other bounds on / and M that 
could come from astrophysics and cosmology. Work is already in progress 
in these directions. In particular, concerning cosmology, it is interesting to 
notice that the allowed range of parameters suggests that branons could have 
weak couplings and large masses. This makes them natural dark matter can- 
didates. In fact, an explicit calculation shows that their relic abundance can 
be cosmologically relevant and could account for the fraction of one third of 
the total energy density of the universe in the form of dark matter presently 
favoured by observations. These results will be presented elsewhere [21]. 
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A Branon vertices 



We show the Feynman rules with leaving momenta for massive branons, in- 
cluding the interaction vertices between two branons and the SM particles. 
The dependence on the momenta of the different particles is explicitly writ- 
ten. We have used the notation in [6]. Different expressions in the case of 
massless branons have been derived in [11]. In that reference the fermionic 
equations of motion are imposed on the Feynman rules. In any case, as a 
consistency check, we have seen that our final results are the same in both 
cases. 



A.l V1[pi,P2,P3,Pa] 



vr"(p3) 



TT' 



\ 



\ 



\ 



\ 



\ 



\ 



/ 





^(Pi) 



VI 



= {7>4a.(P3, Pi - P2) + YPSf^iPi, Pi - P2) 

- nPi,-P2,)ilM^ + 2ips,p,)) 
+ 4m^((p3,P4) + M2)}. 



(54) 
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A.2 V2li[pi,p2,P3,P4] 




+ PliyP4l^{P2,P3) +P4uP2,,{PuP3) 

- Vt^u{{PuP4){P2,P3) + {Pl,P3){P2,P4) " (Pl , P2) (Ps, P4) ) 

- iPl,P2)iP4uP3t^ + P3uPif,) - (P3,P4)(P1/.P2m) 

- ^^aC^P^i^P^i^ + 2p3i.P4/. - '2Vtiu{P3, Pa) - rif,uM^)}, 

(55) 

where we have used the flat background metric to contract indices in the 
Fadeev-Popov Lagrangian. 
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A.3 V31[ps,Pa] 




1p{p2) 



^3)1 = {2YP4uP3,^ + 2YP3uP4,^ 

+ 7^(-3M2-4(p3,p4)}(cy-CA75)- 



(56) 



A. 4 Vr^l\[pi,p2,P3,P4,P5 



7r"(p3) 



vr^(p4) 




yKi\[PuP2,P3,P4,P!i] 



AI{P2) 



Al{p,) 



{Plu {P3^lP4X + PZXPi^i) - PlX {PSfiPiu + P3uP4fi) (57) 



+ P2x{P3uP4,i + P3fj.P4u ) - P2i,{P3uP4X + P3XP4v) 

+ P5m {P3uP4X + P3XP4v ) " Pbu {P3txP4X + P3AP4m ) 

+ ^lXv{{P3,P4){p2^Ji - Pb^l) + P4i,{p^ - P2,P3) + P3^,{P^ - P2,P4)) 

+ ^Am((P3,P4)(P5i' -Pli.) +P4u{,Pl -Pb^Pa) +P3u{Pl -P5,P4)) 

+ Vliv{.{.P3,P4){.PlX -P2X) +P4xiP2-Pl,P3) + P3xiP2 - Pi, Ps))} 
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A. 5 V5[pi,p2,P3,P4\ 



7r°(p3) 



vr^(p4) 



V"5[pi,P2,P3,P4] 



0(P2) 



V5 



7^ 



{-[(P3,P4) + M2][(pi,p2)+m^ 



+ (Pl,P3)(P4,P2) + (P2,P3)(P4,Pl) + -M2(pi,P2) 



30 



A.6 V61[ps,Pa] 



0(^2) 



+ (Pl -P2,P3)P4/. + {Pl -P2,Pi)p^^,]■ 



(59) 



A.7 V7f,[p^,n] 



7r"(p3) 



ah 



—ih'^{T°',T''}6"^ 1 
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A.8 V8li%[ps,p,] 



7r-(p3) 



TT 



■^(P4) 




[P3,P4] 



V8;%i[ps,p4 = '-^{C^'''C^"'[VuAP3pP^,+P3,P4p) 



- Vup{P3aPif, + P3^^Pia) 
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+ VpaiVup{P3, Pa) - P3pPAu - PSuPAp) 
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- VuaiPSpPAp + P3f,PAp) 

+ Vpp iVua {P3 ,Pa) - P3<tPAv " P3uPAa ) 

- VapiVupiP3,PA) - P3p,PAu " P3uPAp)]} ■ 
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B Probability amplitudes 



The probability amplitudes are easily calculated from the Feynman rules 
introduced before: 

B.l Ml: ilj{pi), i/j{p2) 7r(A:i)7r(A:2)^«(g) 

There are four different diagrams that contribute to the tree-level amplitude 
of this process: 




Vj/ Tl \|/ V|/ \|/ 



which is given by the following expression: 

Ml\p^,p2,h,k2,q] = v[pl]{V3;[h,k2]e^'-[q,al] (62) 
+\/li;:£^''^[g, al]tA^j)[p, - q]Vl[-ip, - q), -p2, h, ^2] 
+Vl[~p^, -{p2 - q), h, k2]tA^f)[~{p2 - q)]Vir^e^'^[q, al] 
+\/ir,^Aj;;f [(pi + P2)]V2f^[q, -iPi + P2), h, k2]e^'^[q, cTl]}u[p2]. 

The fermion (iA(j)) and gauge boson (iA(A)) propagators have the usual 
expression in the Standard Model in renormalizable gauges, and we have 
used the standard vertex between gauge and fermions fields: 
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B.2 M2: Al{q) 



If there are massive gauge bosons in the theory, branons contribute to their 
decay through these processes, whose probabihty amphtude can be related 
to the previous one by: 



where in the right hand side, the spinors have been changed as follows ^Ipi] —>■ 
u\pi\ and u\p2\ v\p2\. 

C Integrals over two-body phase space 

Integrals over the two-body phase space of branons can be easily performed 
using: 



M2[pi,p2, h, h, q] = Ml[-pi, -p2, h, k2, -q], 



(63) 




k2)klfjkii,k2pk2a) 

(64) 



+ 



+ 



Cik^kiikpka + C2k'^{rjij,i,rjpa- + rjuprjiia + Tip,aVvp) 

C^k {kij^ki,T]p(^ + kpkcrT]^i,) 

C^k {kjjkprjucr + kiikcrfj^ip + kfjikaflpv ~l~ kpki,rjncr) 



I. 




C^kpki, + C^k'^rjui, 



(65) 




(66) 



And the same result by changing ki for /c2- 




Co 



(67) 



where 
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C^2 = rT\/l TT-ffc -4M^)^ 



D Polarized cross sections 

In this section we give the expressions for the cross sections of the process: 

a A : VS(Pi),V'(P2) n{h),n{k2),Al{q) (69) 

for the different polarizations of the outgoing gauge field. We have used the 
definition: 

— (q'o, |"^| sin^^cos0, sin6'sin0, cos6'), (70) 
The polarization vectors read: 

Sq = —(I'g^l, go sin 6* cos 0, go sin 6* sin 0, go cos 6'); (71) 
e± = ^=(0, — cos^cos^ ± isin0, — cos^sin0 ^ icos0,sin^). 
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For the longitudinal polarization (eg), we get: 
daA _ \h\' {cl + c\)Mln 



^ 4M2 



dxd cose An 122880/%2(s - M^)^\\ s{l - x) + 

yjs{sx^ - AM^){5M\2{s{l - x) + Ml) - hM'^){sx'^ shi^ 9} 

(5(1 - x) + Ml - AM^y 

^ (4M2(M2 + s(l - x)) + s(sx^ - AMI) sin^ 

{32(s(l -x) + Mlf{Ml + sf - 16(s(l -x) + Mlf 

[Ml + 2Mlsx{l - x) + M^s2(9 _ 2^2) 
+ 2s^(l -X- x'^)] siv? - As[AMl + s(12 - 20x + 15x^) 
- My{4 + 8x-37x'^ + llx^ + 4x^) 
+ Mls^{-12 + 2Sx - 7x^ - Gx^ - lOx^ + Ax^) 
+ /a;^(3 - 7a; + 2a;^ + 2x^)] sin^ ^ 

+ 2s\2 + x^){sx^ -AMlfsiTL^e}}. (72) 

The contributions of the transverse polarizations (e^ or e~) are the same and 
they are given by: 



da A {cl, + c\)n 



AM^ 

1 - 



dxdcosO 47r 122880/%2(s_M2)2'^ s{l - x) + Ml 

^s{sx^ - AMl){10M^M^{2{s{l - x) + Ml) - 5M^){2 - sin^ 6} 
(s(l -x) + Ml - AM'^)'^ 
^ {AMl{Ml + s{l - x)) + s{sx^ - AMI) sin^ 

{16M^(s(l - x) + )2[9M^ + s(6 - Ax) + 8^(5 - 4x + 2x^)] 
+ 2[-36Ml^ + 2Ml^s(-108 + A0x + x^) 

- AMls^(68- 72x + 6x^ + x^) 

+ Mls^(-136 + 200x - 54x^ - 26x^ + 3x^) 

- My(AA - 80x + 3Ax^ + 12x^ -13x^ + x^) 
+ Mls^x(-8 + 28x - 3Ax^ + 21x^ - 6x^) 

- s^x^{-2 + 4x - 3x^ + x^)] sin^ 9 + s[120M„^° 

- 2Mf s(-148 + 76x + I3x^) - Mls^(-8 + A8x + 62x^ 

- 34x3 + x^) + M^s=^(24 + 8x + lOx^ + Ax^ - 3x^ + x^) 
+ M2sV(-18 + lOx - 11x2 ^ 2x^) + s^x^(3(l - x) 

+2x2)] sin^ e - s'^isx'^ - AMI)^[AMI - 2Mls + s^x^] sin^ 9}} (73) 
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